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Abstract 

We illustrate how the mathematical modelling of the equations of motion 
in terms of autonomous Hamiltonian dynamical system can definitively fix a 
sign for an otherwise indefinite sign of a certain astrophysical parameter. To 
illustrate it, we shall consider the Mannhcim-Kazanas-de Sitter solution of Weyl 
gravity containing the parameter 7, which is believed to be significant in the 
halo gravity. The strategy we adopt is to calculate the maximum radius up to 
which the halo supports stable material circular orbits. The maximum radius for 
several observed lenses are calculated for both signs of 7, and with the observed 
value of cosmological constant A. These lenses (all having approximately the 
Einstein radius i?E ~ 10 23 cm) consistently yield a maximum radius fl^ax (— 
4.25 x 10 27 cm) inside the de Sitter radius of the universe only when 7 is negative, 
while a positive 7 yields Rn^° always exceeding the de Sitter radius. 
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I. Introduction 

The issue of dark matter (halo matter), arising out of reconciling known 
gravitational laws with observed fiat rotation curves, is a challenging problem in 
modern astrophysics. Observationally, so far, mass of the dark matter appears to 
increase with distance in galaxies, but in clusters exactly the reverse is true, the 
dark matter distribution actually decreases with distance. Indeed, for certain 
dwarfs (such as DD0154) the rotation curve has been measured to almost 15 
optical length scales indicating that the dark matter surrounding this object 
is extremely spread out (see, for a review, Sahni [1]). The total mass of an 
individual galaxy is still somewhat of an unknown quantity since a turn around 
to the v cc r~ x / 2 law at large radii has not been convincingly observed. On 
the other hand, dark matter is attractive and localized on scales less than the 
cosmological distances where repulsive dark energy prevails. Therefore, one 
would expect that maximum radii of stable circular material orbits in the halo 
be smaller than the de Sitter radius of the Universe. Imposing this condition on 
a certain model theory with an undetermined parameter 7, we obtain its sign 
unambiguously. This is the idea in this paper. 

Although classical Einstein's general relativity theory has been nicely con- 
firmed within the weak field regime of solar gravity and binary pulsars, observa- 
tions of flat rotation curves in the galactic halo still lack a universally accepted 
satisfactory explanation. The most widely accepted explanation hypothesizes 
that almost every galaxy hosts a large amount of nonluminous matter, the so 
called gravitational dark matter, consisting of unknown particles not included 
in the particle standard model, forming a halo around the galaxy. One of the 
possibility could be that these particles (WIMPs) encircling the galactic center 
are localized in a thick shell providing the needed gravitational field [2]. The 
exact nature of either the dark matter or dark energy is yet far too unknown be- 
yond such speculations. There also exist alternative theories, such as Modified 
Newtonian Dynamics (MOND) [3,4], braneworld model [5], scalar field model 
[6] etc that attempt to explain dark matter without hypothesizing them. A 
prominent model theory is Weyl conformal gravity and a particular solution in 
the theory is the Mannhcim-Kazanas-de Sitter (MKdS) metric [7] that we shall 
consider here. 

The MKdS solution contains two arbitrary parameters 7 and k (— A/3) that 
are expected to play prominent roles on the galactic halo and cosmological scales 
respectively. While the value of the cosmological constant A = 1.29 x 10~ 56 
cm~ 2 is well accepted, there is some ambiguity about the sign and magnitude 
of 7. From the flat rotation curve data, Mannheim and Kazanas fix it to be 
positive and being of the order of the inverse Hubble length, while Pireaux [8] 
argues for I7I <~ 1CP 33 cm -1 . Edery and Paranjapc [9] fix it to be negative from 
the gravitational time delay by galactic clusters while the magnitude is still of 
the order of inverse Hubble length. Thus there exist great ambiguity both in 
magnitude and sign of 7. 

In the present context, we recall the observational fact that massive neutral 
hydrogen atoms are executing circular motions in the halo around the galactic 
center [10]. The redshiftcd light from those atoms are measured to determine 
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their tangential velocities [11]. Therefore it is necessary to consider massive test 
particle orbits. However, due to conformal invariance of the theory, geodesies 
for massive particles would in general depend on the conformal factor il 2 (x) 7 
but here we assume a fixed conformal frame and not considering other conformal 
variants of the metric. 

The purpose of this article is to illustrate how application of a Hamilto- 
nian system can yield the maximum radius around a galactic center within 
which there can exist stable circular orbits of massive test particles and beyond 
which the orbits become unstable. The criterion of stability has been originally 
suggested by Edery and Paranjape [9] because it provides a way to the determi- 
nation of a natural length scale or region of influence of localized sources in the 
cosmological setting. The strategy we adopt here is to frame the geodesic equa- 
tion in the MKdS solution as a Hamiltonian system, and based on it, analyze 
the stability of circular motion, which then favors a negative 7. 

II. Geodesic equations 

The MKdS metric is given by [7,9] (in units G = 1, vacuum speed of light 
co = 1): 

dr 2 = B(r)dt 2 1 —dr 2 -r 2 (d9 2 +sin 2 9d<p 2 ), B(r) = 1- — + 7 r-fcr 2 , (1) 

B(r) r 

where k and 7 are constants. 

Using u = 1/r, we get the following path equation for a test particle of mass 
mo on the equatorial plane 9 = tt/2 (see Appendix A): 

d 2 u , 7 M If 2k\ 



u + 'iMu 2 -l. (3) 



dv 2= - u + 3Mu -2 + 7^ + <w^--J' (2) 

where h = the angular momentum per unit test mass. For photon, m a = 
=> h — >• 00 and one ends up with the conformally invariant equation but 
without k making its appearance: 

d 2 u „ „ , 9 7 

— = -u + 3Mu 2 - 1 
d^p 1 2 

In the Schwarzschild-de Sitter (SdS) metric, such a disappearance has been 
noted for long [12] but here we find that it occurs despite the presence of 7 
in the metric. The impact of k and 7 on light bending has been investigated 
elsewhere, in Rcfs. [13-16]. 

Analysis of dynamical system involves converting the second order equation 
into two first order equations [17,18] (see Appendix B). For this purpose, we 
introduce the notation 

dx r*\ 
u = x, y = x = — (4) 

to reduce Eq.(2) into a pair of first order autonomous system in the (2, y) phase 
plane 

x = X(x,y)=y (5) 
y = Y(x, y) = a + bx + cx 2 + dx~ 2 + ex~ 3 (6) 
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where 

M t t k 

a= 1? -lb=-l t c = W t d=+ i ,e = - w (7) 

First we discuss stability of circular orbits of light, though their stability is 
not essentially needed. Recall that even in the Schwarzschild spacetime, circular 
light orbits at R — 3M are unstable. Nonetheless, it is instructive to have a 
look at this aspect in the MKdS solution. 

III. Dynamical autonomous system: Massless particle motion 
Light motion occurs in circular orbits defined by R(2 + ^R) — 6M = [see 
Eq.(ll) below] because h 2 — > oo, which implies that d = e = but 7^0. 
The equilibrium points are given by x — 0, y — 0, which yield ^ -NVfr'^ ac ^ Q^j 
and ^ -b-Vb'^^ 4ac^ Q^j rp Q j oca ^- e these points on the real phase plane (x,y), we 

must have a 2 = b 2 - Aac = 1 + 67M > 0. Now a 2 = => 7 = — gjy, so the 
equilibrium points reduce to one single point given by P : (^,0). For a 2 > 0, 
or 7 > — giy, there are two distinct equilibrium points points Q± : (73^-, 0) 
where a = + 67M. Thus Q± correspond to two 7— dependent light radii 
R ± = j-j— Q^== ; which expand as follows 



R+ = zl±VHm wM f + Q( 7 ), (8) 
7 

fl _ = -l- v ^T67^ Rj _ 3M _2 + 

7 7 



We have from Eq.(23) below 



.10, 



which yields 90+ = — VI + 6A#7 < leading to unstable radius at R = R + , 
while go- = \/T + 6M7 > showing that i? = i?_ is a stable radius. The 
basic constraint (reality condition) for both is that 7 > — With, say, 7 = 
—7 x 1CT 28 cm -1 (value inspired by Ref.[9]), the constraint is always satisfied 
for known lenses (say, for M — 2.9 x 10 18 cm, Abell 2744). From Eq.(9), then 
we get the value R- = 2.86 x 10 27 cm, at which there is stability. However, 
stability of massive particle circular orbits is more relevant, which we examine 
next. 

IV. Hamiltonian system: Massive particle motion 

The equilibrium points are given by x — and y = 0. The equation x = 
gives r = R = constant, while y = gives 

2MR 2 + RH 1 -2kR) 

R(2 + 7 R)~6M ■ 1 ' 

The autonomous system (5), (6) can be phrased as a Hamiltonian system as 
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follows 

dU 

— = -Y (x, y) = -(a + bx + ex 2 + dx~ 2 + ex~ 3 ) (12) 
r) H 

— =X(x,y) = y. (13) 

The necessary and sufficient condition for the system (12), (13) to be a Hamilto- 
nian system, namely, ^f- + ^ — 0, is fulfilled for all x and y. Moreover, ^ = 
and therefore H (x, y) — constant (independent of ip). Integrating Eqs.(12),(13), 
we get 

H(x, y) = -{ax + h -x 2 + C -x 3 - dx' 1 - ^x~ 2 ) + u(y) (14) 

H(x,y) = hj 2 + v(x) (15) 

where u(y) and v(x) are arbitrary functions subject to the consistency of Eqs.(12) 
and (13). These two equations will match only if 

u(y) = \y 2 + C (16) 

v(x) = -(ax + ^x 2 + ^x 3 - dx^ 1 - 7^x~ 2 ) + E (17) 

where C, E are arbitrary constants. The family of Hamiltonian paths on the 
phase plane are given by 

H(x, y) = \y 2 - (ax + h -x 2 + °-x 3 - dx^ - ^x~ 2 ) + G (18) 

where G is a parameter. It follows that 
d 2 H 



= -(b + 2cx- 2dx~ 3 - 3ex- 4 ) (19) 

= 1 (20) 
= 0. (21) 



d 2 H 
dy 2 

o 2 h 



dxdy 

As before, the equilibrium points occur when X = and Y = 0, which give the 
values of the orbit radius r = R = constant and h 2 as in Eq.(ll). The quantity 
determining stability is [17] 



_ ( d2H \ 2 



d 2 H d 2 H 

dx 2 dy 2 \ dxdy J 



(22) 



Putting the value of h 2 , we get at the equilibrium points, suffixed by zero, the 
following expression 

6M R(3kR- 7 )[R(2 + 1 R)-6M} 

qo ~ 1 ~^ + i?2 (2 fci?- 7 )-2M • (23) 
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This expression is more general than that in the Schwarzschild spacetime, where 
k = 7 = 0. When q > at any point P : (R,0), we say that P is a stable 
center, but it is an unstable saddle point if q < 0. When q = 0, it is an 
inflection point where the system begins to become unstable. Thus qo = gives 
R = i?max' C 7 beyond which the orbits begin to become unstable. There is also 
a singular radius R — i? s in g where qo blows up. In all cases we have studied, 
-Rsing < -Re for 7 negative and i? s ing > -Rds f° r 7 positive. In either case, the 
occurrence of such singularity is a bit queer but fortunately it occurs at radii 
not accessible to gravitational lensing observations. 

We shall calculate the values of qo for observed lenses. Due to advances 
in technology, mass M and Einstein radius Re for several lenses are available. 
The observed lens masses M are used to compute the values of q marching 
from the Einstein radius Re out to i?ds- A meaningful stable radius should 
lie between Re and i?ds = a/A/3 = 1.52 x 10 28 cm. Fig.l illustrates for one 
of the observed lenses (Abell 2744, M = 2.90 x 10 18 cm, i? E = 2.97 x 10 23 cm), 
which shows that stable material radii exist for all radii R > Re extending 
even beyond the dS radius when 7 = +7 x 1CP 28 cm -1 , which seems physically 
unlikely due to the fact that in the regions beyond dark matter, repulsive dark 
energy takes over. In the cosmological region, where the spacetime becomes 
time dependent and accelerating, the existence of circular orbit itself becomes 
impossible. Also, the singular radius occurs at i? s in g = 8.14 x 10 28 cm > i?ds- 
(Fig. 2). On the other hand, when 7 = —7 x 10~ 28 cm -1 , stability is achieved 
up to a radius R = 4.25 x 10 27 cm, which is well within i?ds and beyond that 
maximum radius instability begins to appear. In this case, the singular radius 
occurs at i? s ing = 9.11 x 10 22 cm < Re (Fig. 3). Thus the maximum radius up 
to which stable material circular orbits are admissible is i?max lc = 4-25 x 10 27 
cm for a negative 7. This limit is not much sensitive either to the exact value 
of A or to other negative values of 7 (such as argued by Pireaux [8]). It can be 
verified that, for other lens data as well (as listed in Ishak et al. [14], but not 
tabulated here), the maximum stable radius remains very nearly the same. We 
must caution that this -R^ax' i s n °t the size of the halo because it would be too 
big; it is just the maximum allowed radius of stable orbits. 

V. Conclusion 

Dynamical systems have been in use in gravitation theories, particularly in 
the regime of cosmology. The application of the dynamical system approach 
to higher order gravity has been shown to be a powerful tool, revealing global 
features of curvature quintessence models (see, for instance, [21]), where the 
cosmological equations per se are posed as a dynamical system in higher di- 
mensional phase space. Posing the geodesic equations as dynamical system for 
studying their stability within the galactic halo regime also seems interesting in 
its own right. 

From the foregoing analysis, we conclude the following: A finite maximum 
stable radius R!^^ c exists within i?ds only when 7 is negative. Much less is yet 
conclusively known about the halo content or size, which has engendered many 
theoretical models. The occurrence of singularity in R is actually a defect of 
the model under consideration, not of the Hamiltonian method. The method is 
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quite general and any other metric solution, if available for the halo, can be used 
equally. The key physical input we used was that the stable material (neutral 
HII) finite circular orbits exist only inside the halo (giving rise to flat rotation 
curves [10]) and not outside it. The Hamiltonian analysis supports such orbits 
only for a negative 7. This is the main result of the paper, which is probably is 
of some interest in view of the fact that it potentially decides the long standing 
ambiguity in the sign of 7 mentioned earlier. 
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Appendix A 
The standard form of geodesic equation is given by 



dx v dx x 
g v \dx v dx x — dr 2 . 



i_ = n 

dr 2 c?t dr ~ ' 



Using the expression for Christoffcl symbols , it can be easily rewritten in a 
convenient form as (see, for instance, [20]): 

dU E _l_d9^ u » u x =0 (A1) 
dp 2 dx» v ' 

g„ x U»U x = m 2 , (A2) 

where dr = m dp, m is the rest mass of the test particle, U x (— ^-) is the 
4-velocity, here x M = (t,r,9,ip) and g v \ is the metric tensor given by (1). The 
advantage in redefining the affine parameter r — > p is that the equations (Al), 
(A2) are now applicable to both massive and massless particle motion since 
dp ^ 0. 

When we construct the (p = 2) 9— component from (Al), we find that 
9 must be a constant, that is, motion takes place in a plane exactly as in 
the Newtonian central force problem. By suitable orientation of axes, we can 
set 9 = n/2. Next, since the metric tensor g u \ do not depend on t and ip, 
the corresponding components of must be the constants of motion. These 
constants are Uq and U3 given by 

U =gooU° = B(r)^, (A3) 
dp 

t/ 3 =.933t/ 3 = -r 2 ^. (A4) 
The remaining constant is just the rest mass mo. Then from (A2), we obtain 

(^) 2 = UZ-U!^- m lB(r). (AS) 
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Using (A3), we can rewrite (A5) as 



*y-BV)-f£V£M-$*(r). (A6) 



dt J w \U J r 2 U$ 

From (A3) and (A4), we get 

dip = U 3 B(r) 
dt U r 2 ' 

Eliminating t between (A6) and (A7), we obtain 
drV [U^ 2 



(A7) 



*p) ={W 3 ) ^-^W-^(r). (A8) 
Putting B(r) = 1 — — + -fr — kr 2 and setting u= 1/r, we can rewrite (A8) as 



du\ 
dip] 



U 2 • \U Q 



- v 2 + 2Mu 3 - U7 + 



2Mm,QU km 2 , jrn 2 , 
U 2 + uHJ 2 ~ uU 2 ' 



The term in square bracket on the right is constant. Differentiating both sides 
of the above with respect to ip and defining ft = we get 

d 2u 2 7 M 1 / 2fc\ /A ^ N 

U + 3M« 2 - J + + 7-- , (A9) 



dip 2 2 h 2 2h 2 u 2 V u 

which is Eq.(2) in the text. 

Appendix B 

It is useful to state the type of (in)stabilities of an equilibrium state we 
had been discussing in the text. An equilibrium state corresponds to a con- 
stant solution of a differential equation describing a mechanical system and 
conversely. Constant solution means that velocity x(= y) and acceleration x 
be simultaneously zero. The concept of (in)stability of an equilibrium state is 
best understood from the familiar example of motion of a pendulum about the 
equilibrium state x — and x = 0, where x is angle with the vertical. The 
motion is stable because a small displacement from the equilibrium state will 
lead to only small oscillations of the bob about the position (0, 0) in a vertical 
plane. 

An equilibrium state in actual space is represented by a point on the phase 
space (x,x). Different closed phase paths y = f(x) around an equilibrium point 
on the phase plane (x, y) correspond to oscillations of the bob in the actual (x, t) 
space with different periods. In the dynamical terminology, the point (0, 0) on 
the phase plane is then called a stable center. If a small displacement from any 
equilibrium state takes the system far away from it, the state is represented by 
an unstable equilibrium point on the phase plane. For instance, the equilibrium 
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state x — 7r and x — 0, that is, the point (tt, 0) on the phase plane is an unstable 
saddle point in the pendulum motion because due to a small perturbation, the 
bob moves away from the point (tt, 0) never to return. This is reflected on the 
phase plane by open paths originating from the equilibrium point. Finally, there 
could be whirling motions of the bob off the vertical plane. This is picturized 
on the phase plane by open paths never passing through an equilibrium point 
and there is no question of stability here. 

With the above understanding, the dynamical framework proceeds with lin- 
earizing the general system [17] 

x = X(x,y) (Bl) 
y = Y(x,y) (B2) 

around an equilibrium point (xq, yo) such that 

X{x ,y ) = 0, (B3) 
Y(x ,y ) = 0. (B4) 

(Note that the family of phase paths y = f(x) is obtained by integrating ^ = 
y)- Taylor expanding around (xo,j/o) retaining only first order derivatives, we 
obtain from (Bl), (B2) the linear system 

x = a(x - x ) + b(y - y ), y = c(x - x ) + d(y - y ), (B5) 

where 



a b 




dx 


dX 




dx 




c d 




dY 




dx 


dy 



the derivatives being calculated at (xo, yo). Depending on the values of p, q and 
A defined by 

p = a + d, q = ad — be, A = p 2 — Aq, (B7) 

one defines the stable center (p = 0,q > 0, A < 0), unstable saddle point (q < 
0, A > 0), stable node (p < 0, q > 0, A > 0), stable spiral (p < 0, q > 0, A < 0), 
unstable node (p > 0, q > 0, A > 0), unstable spiral (p < 0, q > 0, A < 0), 
degenerate stable node (p < 0, q > 0, A = 0), degenerate unstable node (p > 
0,q > 0, A = 0) etc. These classifications, based on the eigen values of the 
matrix (B6), characterize on the phase plane the paths in the neighborhood of 
(xo,yo). 

It is known that the characterization of a stable center via linear approxi- 
mation alone is 'fragile' [17], hence unreliable, due to the fact that it demands 
a rather stringent exact relation (a + d = 0) to be fulfilled by the first order 
derivatives. There are indeed specific problems where a stable center defined 
in this way reflects only a transition between stable spiral and unstable node 
whereas the true behavior of the full system (Bl), (B2) is entirely different. A 
more reliable approach is to consider the extrema of the Hamiltonian H defining 
a center (q > 0), as we did in the text. If we put, from Eqs.(12) and (13), the 
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relations = X (x, y) and ^ = — Y (x,y) both evaluated at (xo,yo), into 
(B7), we immediately see that p — a + d = identically. Also, the value of q at 
(:ro, 2/o ) becomes exactly what we designated in the text as qo- Clearly, this is a 
second derivative approach, hence conclusiv^- When qo = 0, the extrema of H 
marginally corresponds to a center but complicated saddles are not excluded. 
When qo < 0, the equilibrium point (xo, yo) is definitely unstable, likely to be a 
saddle. 

Observationally, we know that there are material orbits in the halo reflecting 
the light coming to us and those orbits must be stable to perturbations. Other- 
wise, such orbits would have collapsed to the galactic center. Thus they should 
correspond to stable center representing periodic motion in the real (x, t) plane, 
which is the physical criterion we used. However, stable nodes or spirals on 
the phase plane are not ruled out in principle, but none of them represents any 
periodic motion in the real (x, t) plane as opposed to what is observed in the 
redshift measurements. The speculation of huge amount of non-luminous dark 
matter being hidden in the halo originates from the observed Doppler emissions 
from stable circular orbits of neutral hydrogen clouds [2,5,10,11,19]. Therefore, 
stable modes on the phase plane other than that of the stable center should be 
ruled out on observational grounds. 

Figure captions: 

Fig.l. The upper line corresponds to 7 = 7 x 10~ 28 cm -1 and the lower line 
to 7 = —7 x 10~ 28 cm -1 in the case of Abell 2744. The plot remains essentially 
the same for other lenses tabulated in Ref. [14]. 

Fig. 2. Curves I correspond to 7 = 7 x 10~ 28 cm -1 giving a singularity of 
qo at R sing = 8.14 x 10 28 cm > R ds . 

Fig. 3. Curves II correspond to 7 = —7 x 10~ 28 cm -1 giving a singularity 
of q at i? S ing = 9.11 x 10 22 cm < R E . 
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